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Designed as a high-sensitivity gamma-ray observatory, the Fermi Large Area Telescope is also an 
electron detector with a large acceptance exceeding 2 m'^sr at 300 GeV. Building on the gamma- 
ray analysis, we have developed an efficient electron detection strategy which provides sufficient 
background rejection for measurement of the steeply-falling electron spectrum up to 1 TeV. Our 
high precision data show that the electron spectrum falls with energy as E^^'" and does not exhibit 
prominent spectral features. Interpretations in terms of a conventional diffusive model as well as a 
potential local extra component are brieffy discussed. 

PACS numbers: 96.50.sb, 95.35.+d, 95.85.Ry 



Introduction. - Accurate measurements of high-energy 
Cosmic Ray (CR) electrons {we hereafter refer to elec- 
trons as a sum of e+ and unless specified otherwise) 
provide a unique opportunity to probe the origin and 
propagation of CRs in the local interstellar medium and 
constrain models of the diffuse gamma-ray emission [1]. 
Prior to 2008, the high-energy electron spectrum was 



'Electronic address: luca.latronico@pi.infn.it 

t Electronic address: AIexander.A.Moiseev@nasa.gov 



measured by balloon-borne experiments [2] and by a sin- 
gle space mission (AMS-01, [-i]). The measured fiuxes 
differ by factors of 2 — 3. While these data allowed sig- 
nificant steps forward in understanding CR origin and 
propagation, constraints on current models remain weak. 
The CR propagation package GALPROP [4], on the as- 
sumption that electrons originate from a distribution of 
distant sources mainly associated with supernova rem- 
nants and pulsars, predicts a featureless spectrum from 
10 GcV up to few hundreds of GcV. Above that energy, 
due to the actual stochastic nature of electron sources in 
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space and time, and to the increasing synchrotron and 
inverse Compton energy losses, the spectral shape may 
exhibit spatial variations on a scale of a few hundred 
parsecs. Nearby sources start contributing significantly 
to the observed local flux and may induce important de- 
viations from a simple power law spectrum [2, 5, 6]. 

Recently published results from Pamela [7], ATIC [8], 
H.E.S.S. [9] and PPB-BETS [10] have opened a new 
phase in the study of high energy CR electrons with a 
new generation of instruments. These four experiments 
report deviations from the reference model mentioned 
above. Pamela measures an increase of positrons with 
respect to electrons at energies above a few GeV; ATIC 
and PPB-BETS detect a prominent spectral feature at 
around 500 GeV in the total electron plus positron spec- 
trum; H.E.S.S. reports significant steepening of the spec- 
trum above 600 GeV. All these results may indicate the 
presence of a nearby primary source of electrons and 
positrons. The natures of possible sources have been 
widely discussed, two classes of which stand out: nearby 
pulsar(s) ( [11], [12] and references therein) and dark mat- 
ter annihilation in the Galactic halo, e.g. [l.'-!]. 

The Large Area Telescope (LAT) is the main instru- 
ment on-board the recently launched Fermi Gamma-Ray 
Space Telescope mission. It is conceived as a multi- 
purpose observatory to survey the variable gamma-ray 
sky between 20 MeV and 300 GeV, including the largely 
unexplored energy window above 10 GeV. It is therefore 
designed as a low aspect ratio, large area pair conversion 
telescope to maximize its field of view and effective area. 
The LAT angular, energy and timing resolution rely on 
modern sohd state detectors and electronics [14]: a multi- 
layer silicon-strip tracker (TKR), interleaved with tung- 
sten converters with a total depth of 1.5 radiation length 
(Xq) on-axis; a hodoscopic CsI(Tl) calorimeter (CAL), 
8.6 Xq deep on-axis; a segmented anticoincidence plastic 
scintillator detector (ACD) and a flexible, programmable 
trigger and filter logic for on-board event filtering. 

Since electromagnetic (EM) cascades are germane to 
both electron and photon interactions in matter, the 
LAT is also by its nature a detector for electrons and 
positrons. Its potential for making systematics-limited 
measurements of CR electrons was recognized during the 
initial phases of the LAT design [l-'j], [[(>]. In this Letter 
we describe the technique developed for this measure- 
ment and our validations using ground and flight data, 
and present the first high statistics CR electron spectrum 
from 20 GeV to 1 TeV based on the data taken in the first 
six months of the mission. A more complete description 
of the analysis procedure, and new results with increased 
statistics, energy coverage and a preliminary study of 
sensitivity to anisotropics in the distribution of arrival 
directions will be presented in forthcoming papers. 

Event selection. - An analysis that separates elec- 
trons from the dominant CR hadrons is required to fully 
exploit LAT's large collecting power and long observa- 



tion time. The on-board filter is configured to accept all 
events that deposit at least 20 GeV in the calorimeter; 
thus we ensure that the rare high-energy events, includ- 
ing electrons, are available for thorough analysis on the 
ground. We developed a dedicated event selection for 
high energy electrons that provides a large geometry fac- 
tor with a residual hadron contamination less than 20% 
at the highest energy (see table I). At any given energy, 
the geometry factor (GF) is defined as the proportional- 
ity factor relating the rate of events passing the selection 
criteria to the incident fiux. It measures the instrument 
acceptance, i.e the integral of the effective area over the 
instrument field of view. It is numerically evaluated using 
a MC data sample of pure electrons. As for the analy- 
sis developed for extracting LAT photon data [14], the 
electron selection essentially relies on the LAT capabil- 
ity to discriminate EM and hadronic showers based on 
their longitudinal and lateral development, as measured 
by both the TKR and CAL detectors. The background 
rejection power for photon science is optimized up to 300 
GeV. The electron selection criteria are instead tuned in 
the multi-100 GeV range, where the much steeper elec- 
tron spectrum requires an overall hadron rejection power 
of 1 : lO''. 

Events considered for the electron analysis are re- 
quired to fail the ACD vetoes developed to select pho- 
ton events [ i]. This removes the vast majority of the 
potential gamma-ray contamination. In fact, the geome- 
try factor for photons, determined with electron selection 
cuts, is less than 8% of that for electrons at 1 TeV. The 
ratio of photon to electron fiuxes is negligible at low ener- 
gies and rises to around 20% at the highest energy. This 
estimate is obtained from a simple, conservative extrap- 
olation to high energies of the EGRET all-sky spectrum 
at GeV energies using a i?^^-^" power law [17]. The over- 
all gamma contamination in the final electron sample is 
therefore always less than 2%. 

EM showers start developing in the LAT TKR, while 
most of the energy is absorbed in the CAL. The mea- 
surement of the lateral shower development is a powerful 
discriminator between more compact EM showers and 
wider hadronic showers. We select on variables that map 
the distribution of TKR clusters around the main track, 
and in the CAL the truncated, second-order moments of 
the energy distribution around the shower axis. A further 
selection derives from the different distributions of energy 
and hits in the ACD between EM and hadron-initiated 
showers. At this stage, the hadron rejection power is at 
the level of 1 to a few 10^, improving to greater than 10'^ 
below 100 GeV thanks to the ACD selection. 

Similarly to the LAT photon background rejection 
analysis, the remaining necessary boost in the rejection 
power is obtained by combining two probability variables 
that result from training classification trees (CT) to dis- 
tinguish between EM and hadron events [14]. This is 
done using large sets of Monte-Carlo (MC) events gencr- 



4 



ated by the accurate LAT simulation package [14], based 
on the Geajit4 toolkit [is]. Two CTs are used, one built 
with TKR variables, and a second one based on CAL 
variables, which describe the complete event topology. 
The variables given most weight by the CTs are the 
same or equivalent to those described above. The clas- 
sifiers allow selection of the electrons through a multi- 
tude of parallel paths, each with different selections, that 
map the many different topologies of the signal events 
into a single, continuous probability variable that is used 
to simultaneously handle all valid selections. The TKR 
and CAL electron probabilities are finally combined to 
create an energy-dependent selection that identifies elec- 
trons with greater efficiency and optimized background 
rejection with respect to a single sequence of cuts. The 
resulting rejection power is flat and better than 1 : 10'^ up 
to 200 GeV and from there rises steadily to 1 : lO'* at 
1 TeV in a manner that partially compensates for the in- 
creasingly larger relative proton fluxes with energy. Con- 
versely, the electron selection efficiency, calculated as the 
ratio of selected versus triggered events, has a peak value 
of 50% at 20 GeV and steadily decreases down to 12.5% 
at 1 TeV. 

Energy reconstruction and validation. - Energy recon- 
struction is the other critical aspect of this analysis. For 
EM cascades of several hundreds of GeV a large fraction 
of the energy falls outside of the LAT CAL. The shower 
imaging capability is therefore crucial in fitting the lon- 
gitudinal shower profile in order to correct for the energy 
leakage and estimate the incoming energy with good ac- 
curacy. The resulting energy resolution for events passing 
the electron selection is shown in figure 1. Since showers 
are not fully contained above 20 GeV, the distribution of 
the reconstructed energy after leakage correction is asym- 
metric, with a longer tail toward lower energies. For this 
reason we quote the full width of the 68% containment 
of the distribution as our energy resolution, and check 
that the full 95% containment does not imply indefinitely 
long tails; see figure 1. Candidate electrons traverse on 
average 12.5 radiation lengths, resulting from the total 
thickness of the TKR and CAL detectors and the effect 
of event selection. 

The energy reconstruction algorithm and the event 
analysis rely heavily on the LAT MC simulation. This 
was extensively verified and fine-tuned using beam test 
data for electrons and hadrons up to 282 GeV [19]. 
Extensive efforts are made to avoid bias in the event 
selection by systematically comparing fiight data and 
MC distributions of likely discriminants of electrons and 
hadrons, and choosing only those that indicate a good 
agreement. Figure 2 shows the very good data-MC 
agreement for the critical variable that maps the trans- 
verse shower size. 

Systematic uncertainties are determined for all energy 
bins. For each step in the event selection, we scan a range 
of thresholds around the reference value used by the cut 
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FIG. 1: (color online) Energy resolution for the LAT after 
electron selection; the full widths of the smallest energy win- 
dow containing the 68% and the 95% of the energy dispersion 
distribution are shown. The comparison with beam test data 
up to 282 GeV and for on-axis and at 60° incidence shown 
in the figure indicates good agreement with the resolution 
estimated from the simulation. 
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FIG. 2: (color online) Distribution of the transverse sizes of 
the showers (above 150 GeV) in the CAL at an intermediate 
stage of the selection, where a large contamination from pro- 
tons is still visible. Flight data (black points) and MC (gray 
solid line) show very good agreement; the underlying distri- 
butions of electron and hadron samples are visible in the left 
(red) and the right (blue) peaks respectively. 



and derive the corresponding flux versus GF curve. We 
extrapolate the curve to a GF consistent with a null cut, 
and take the relative difference of the corresponding flux 
and the reference as the systematic uncertainty associ- 
ated with the cut. All such contributions, taken sepa- 
rately with their signs, and the uncertainty of the resid- 
ual contamination, derived from an overall 20% uncer- 
tainty in the underlying proton spectrum are summed in 
quadrature. The result is shown in table I. 
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The final tuning of the event selection provides a maxi- 
mum systematic error less than 20% at 1 TeV. The abso- 
lute LAT energy scale, at this early stage of the mission, 
is determined with an uncertainty of This esti- 

mate is being further constrained using flight and beam 
test data. The associated systematic error is not folded 
into those above as it is a single scaling factor over the 
whole energy range. Its main effect is to rigidly shift the 
spectrum by ^20% without introducing significant defor- 
mations. 

While event selection is explicitly energy-dependent to 
suppress the larger high-energy background, it is not op- 
timized versus the incident angle of incoming particles. 
Nonetheless we have compared the spectra from selected 
restricted angular bins with the final spectrum reported 
here; they are consistent within systematic uncertainties. 
A further validation of the event selection comes from 
an independent analysis, developed for lower-energy elec- 
trons, which produces the same results when extended up 
to the the endpoint of its validity at ~ 100 GeV. Our ca- 
pability to reconstruct spectral features was tested using 
the LAT simulation and the energy response from fig- 
ure 1. We superimposed a Gaussian line signal, centered 
at 450 ± 50 GcV rnis, on a power law spectrum with an 
index of 3.3. This line contains a number of excess counts 
as from the ATIC paper [y], rescaled with the LAT GF. 
We verified that this analysis easily detects this feature 
with high significance (the full width of the 68% contain- 
ment energy resolution of the LAT at 450 GeV is 18%). 

Results and discussion. - More than 4M electron 
events above 20 GeV were selected in survey (sky scan- 
ning) mode from 4 August 2008 to 31 January 2009. En- 
ergy bins were chosen to be the full width of the 68% 
containment of the energy dispersion, evaluated at the 
bin center. The residual hadronic background was es- 
timated from the average rate of hadrons that survive 
electron selection in the simulations, and subtracted from 
the measured rate of candidate electrons. The result is 
corrected for finite energy redistribution with an unfold- 
ing analysis [20] and converted into a flux Je by scaling 
with the GF, see table I. The distribution of x Je is 
shown in table I and in figure 3. 

Fermi data points visually indicate a suggestive devi- 
ation from a fiat spectrum. However, if we conserva- 
tively add point-to-point systematic errors from table I 
in quadrature with statistical errors, our data are well 
fit by a simple normalized £;~3.04 power law (x^ = 9.7, 
d.o.f. 24). 

For comparison, we show a conventional model [1] for 
the electron spectrum, which is also being used as a ref- 
erence in a related Fermi-LAT paper [21] on the Galactic 
diffuse gamma-ray emission. This uses the GALPROP 
code [4], with propagation parameters adjusted to fit a 
variety of pre- Fermi GR data, including electrons. This 
model has an electron injection spectral index of 2.54 
above 4 GeV, a diffusion coefficient varying with energy 
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FIG. 3; (color) The Fermi LAT CR electron spectrum (red 
filled circles). Systematic errors are shown by the gray band. 
The two-headed arrow in the top-right corner of the figure 
gives size and direction of the rigid shift of the spectrum im- 
phed by a shift of t.^'ox of the absolute energy, corresponding 
to the present estimate of the uncertainty of the LAT energy 
scale. Other high-energy measurements and a conventional 
diffusive model [i] are shown. 



as E^^^, and includes a diffusive reacceleration term. As 
can be clearly seen from the blue dashed line in figure 3, 
this model produces too steep a spectrum after prop- 
agation to be compatible with the Fermi measurement 
reported here. 

The observation that the spectrum is much harder than 
the conventional one may be explained by assuming a 
harder electron spectrum at the source, which is not 
excluded by other measurements. However, the signif- 
icant flattening of the LAT data above the model pre- 
dictions for E > 70 GcV may also suggest the pres- 
ence of one or more local sources of high energy CR 
electrons. We found that the LAT spectrum can be 
nicely flt by adding an additional component of pri- 
mary electrons and positrons, with injection spectrum 
.-/cxtra(£') oc E''''' cxp{~ E / E^nt} , Ecut bciug thc cut- 
off energy of the source spectrum. Thc main purpose 
of adding such a component is to reconcile theoretical 
predictions with both thc Fermi electron data and thc 
Pamela data [7] showing an increase in thc e+/(e~ +e+) 
fraction above 10 GcV. Thc latter cannot be produced 
by secondary positrons coming from interaction of the 
Galactic CR with the ISM. Such an additional compo- 
nent also provides a natural explanation of the steepen- 
ing of the spectrum above 1 TeV indicated by H.E.S.S. 
data [!J]. As discussed in [12] and references therein, pul- 
sars are the most natural candidates for such sources. 
Other astrophysical interpretations (e.g. [22]), or dark 
matter scenarios, can not be excluded at the present 
stage. 

A detailed discussion of theoretical models lies out- 
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146.5 ± 0.71^1 


53.7-60.4 


2.77 


0.12 


194062 


145.5 ± 0.7^71 


60.4-68.2 


2.76 


0.13 


155585 


143.2 ± 0.7tg g 


68.2-77.4 


2.73 


0.14 


126485 


141.9 ± 0.8^7 


77.4-88.1 


2.71 


0.14 


100663 


140.8 ± O.stj'l 


88.1-101 


2.68 


0.15 


77713 


139.0 ± O.gtg g 


101-116 


2.64 


0.16 


61976 


139.0 ± 0.9t^'2 


116-133 


2.58 


0.17 


46865 


139.4 ± l.ot^-j 


133-154 


2.52 


0.17 


35105 


139.5 ± 1.2^7 4 


154-180 


2.44 


0.17 


27293 


140.8 ± 1.3^74 


180-210 


2.36 


0.18 


19722 


142.3 ± 1.5^7:4 


210-246 


2.27 


0.18 


13919 


140.9 ± 1.7j:e'8 


246-291 


2.14 


0.18 


10019 


140.9 ± 1.9lg;7 


291-346 


2.04 


0.18 


7207 


140.4 ± 2.2^7'q 


346-415 


1.88 


0.18 


4843 


139.4 ± 2.6^7 2 


415-503 


1.73 


0.19 


3036 


134.0 ± 3.lt? j 


503-615 


1.54 


0.20 


1839 


127.4 ± 4.1^8 6 


615-772 


1.26 


0.21 


1039 


115.8 ± 4.8+Jo g 


772-1000 


0.88 


0.21 


544 


114.4 ± 6.5tJ?-j 



TABLE I: Geometry factor, residual contamination, number 
of counts before background subtraction and the flux Je mul- 
tiplied by E^. The statistical error is followed by the system- 
atic error. The latter does not include the effect due to the 
uncertainty in the absolute energy scale. 

side the scope of this work and will be presented in a 
forthcoming Fermi interpretation paper. At this stage 
it suffices to say that Fermi will significantly change the 
understanding of this part of the electron spectrum. 
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